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First-principles calculations are performed for electronic structures of two high pressure phases of 
solid selenium, /3-Po and bcc. Our calculation reproduces well the pressure-induced phase transition 
from /3-Po to bcc observed in selenium. The calculated transition pressure is 30 GPa lower than 
the observed one, but the calculated pressure dependence of the lattice parameters agrees fairly 
well with the observations in a wide range of pressure. We estimate the superconducting transition 
temperature Tc of both the /3-Po and the bcc phases by calculating the phonon dispersion and the 
electron-phonon interaction on the basis of density-functional perturbation theory. The calculated 
Tc shows a characteristic pressure dependence, i.e. it is rather pressure independent in the /3-Po 
phase, shows a discontinuous jump at the transition from /3-Po to bcc, and then decreases rapidly 
with increasing pressure in the bcc phase. 

61.50.Ks, 71.15.Mb, 74.25.Jb, 74.25. Kc 



I. INTRODUCTION 

Recently, with the development of high-pressure exper- 
imental techniques, research on pressure-induced struc- 
tural phase transitions of the group VIb elements O, S, 
Se and Te has progressed greatly. At high-pressures both 
metallic tellurium (Te) and selenium (Se) transform from 
/3-Po to bcc phase at 27 GPatl and 150 GPa,H, respec- 
tively. The /3-Po type structure is rhombohedral and can 
be described as a simple cubic lattice deformed along 
the [111] direction keeping the edge length unchanged. 
By changing the ratio c/a of the rhombohedral lattice 
we obtain the bcc structure when c/a — V6/4. In Te 
superconducting transitions are observed in both /3-Po 
and bcc structures at low temperature and a junap in Tc 
from 2.5 K to 7.4 K is observed at 32-35 GPa.l3 Theo- 
retically, Mauri, et. al. suggests that the jump in Tc is 
related to the phonon softening in the bcc phase, namely 
with decreasing pressure the jphonon anomaly enhances 
the electron-phonon coupling.0 In Se, on the other hand, 
there is neither experimental observation nor ah initio 
calculation for the pressure dependence of phonon fre- 
quencies, the electron-phonon interaction or Tc. 

The purpose of the present paper is to estimate these 
quantities by using ab initio calculations. First we calcu- 
late the total energies of the /3-Po and bcc phases of Se 
by using the full-potential linearized muffin-tin orbital 
(FPLMTO) methoda in order to discuss the pressure- 
induced phase transition from /3-Po to bcc. Then we 
calculate phonon frequencies and electron-phonon cou- 
pling constants usins.|the linear-response FPLMTO (LR- 
FPLMTO) method.la The validity of this method is 



demonstrated in references 5-8. Finally we calculate the 
superconducting transition temperature in both the /3-Po 
and bcc phases as a function of pressure. 

II. STRUCTRAL TRANSITION FROM /3-PO TO 
BCC IN SELENIUM 

A. Calculational procedure 

The calculations of electronic states for the /3-Po and 
bcc structures of Se have been done according to the fol- 
lowing procedure with use of the FPLMTO program. For 
exchange-correlation functional we have adopted the for- 
mula proposed by Gunnarsson and Lundqvism and the 
GGA correction proposed by Perdew et. a/JlJ has been 
taken into account. Inside the mufhn-tin (MT) spheres 
the scalar-relativistic calculations are performed for va- 
lence electrons, and the core states are recalculated at 
each self-consistent iteration with relativistic effects. The 
MT radius has been taken to be 1.07 A. The fc-space in- 
tegration hasjpeen performed by the improved tetrahe- 
dron methodliJ with use of (12, 12, 12) grid of the sam- 
pling fc-points [189 points in the IBZ ]. We have used 3k- 
spd-LWIO basis set (27 orbitals): = -0.1 ,-1.0 and 
—2.0 Ryd. In the interstitial region the basis functions 
are expanded in plane waves up to the cutoff correspond- 
ing approximately to 200, 350 and 650 plane waves per s, 
p, and d orbitals, respectively. The charge densities and 
the potentials are expanded inside the MT spheres by 
spherical harmonics up to ^max=6 and in the interstitial 
region by plane waves with the cutoff corresponding to 
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the (16, 16, 16) fast-Fourier-transform(FFT) grid in the 
unit cell of direct space. The final convergence is within 
10-6 Ryd. 



B. Results 

As the first step to investigate the phase transition 
from /3-Po to bee we have calculated the electronic band 
structure and the total energy of f3-Po Se at atomic vol- 
ume Va=14.82A. At this volume the /3-Po structure is 
known to be stable by experimental measurements. For 
comparison we have calculated the electronic band struc- 
ture and the total energy of the hypothetical bcc struc- 
ture with the same atomic volume. The calculated energy 
band structures and the DOS for the /3-Po and the bcc 
structures are shown in Fig. |^. The lowest band is mainly 
derived from the 4s component and the next three bands 
from the 4p component. The band structures of both 
the structures are similar to each other on the whole. 
However, remarkable differences can be seen in the band 
structure and the DOS near the Fermi level. Firstly, the 
band structure of the bcc structure along the P and A 
lines has much larger dispersion than that of the /3-Po 
type structure. Secondly, the DOS at the Fermi level of 
the bcc structure is larger than that of the /3-Po struc- 
ture because in bcc Se the L point is a saddle point of 
the third energy band from the bottom which crosses the 
Fermi level near the L point. By deforming from bcc to 
/3-Po the L point energy of the third band goes up away 
from the Fermi level and the states of the fourth band 
in the middle of the P line come down under the Fermi 
level. As a result the DOS of (3-Po Se has a relatively 
large peak at 0.225 Ryd. below the Fermi level. This is 
the reason why the /3-Po type structure is relatively more 
stable compared with the bcc structure. 

To investigate the pressure-induced structural transi- 
tion we have to calculate the total energy of the (3-Po and 
the bcc structures as a function of volume. For the /3-Po 
structure we have optimized c/a at each volume, namely, 
we have calculated the total energy of the (3-Po type Se 
as a function of c/a with the atomic volume Va being 
kept constant. Figure ^ shows the total energy vs c/a 
for several fixed atomic volume and the inset shows the 
results for a volume range near the phase transition. All 
the energies are referenced to that of the bcc structure, 
i.e. the energy at c/a=\/6/4. For a large volume such 
as 16.30 the energy takes the minimum at c/a ~ 0.85 
and the c/a of bcc is an inflection point. With decreasing 
volume the value of the c/a starts to decrease gradually 
toward the bcc structure. Finally, the (3-Po structure is 
no longer a quasistable structure, that is, the total en- 
ergy has a single minimum at c/ a—\f%/ \ corresponding 
to the bcc structure. 

Figure || shows the total energy of the (3-Po and the 
bcc structures of Se as a function of volume. At larger 
volumes the (3-Po structure is more stable than the bcc 



structure. At smaller volumes the total energies of both 
the structures are quite close, but the bcc structure is 
more stable than the (3-Po structure at volumes smaller 
than ~ 12 A^. 

In order to estimate the transition pressure from (3-Po 
to bcc we calculate the Gibbs free energy (or enthalpy) 
as a function of pressure. Then, to evaluate the pressure 
as a function of volume we fitted the calculated total 
energies by the Murnaghan's equation of state (E0S)li3: 
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where Bq and B'q is the isothermal bulk modulus at zero 
pressure and it's derivative, respectively. The pressure is 
determined from 



B'o 



Vo 



The Gibbs free energy is defined by G(P) = Etot{P) + 
PV{P) and the transition pressure between the two 
phases is obtained from the relation Gp{P) = G'b(P), 
where Gp and Gb are the Gibbs free energies of the (3-Po 
and the bcc structures, respectively. 

The transition pressure Pc from (3-Po to bcc has been 
estimated as 120 GPa by our present calculation. This 
value iSphigher than other. calculated transition pressures, 
90 GPaliil and 110 GPaM^jut still lower than the experi- 
mental value of 150 GPa.ll3 The origin of this discrepancy 
between theory and experiment may be ascribed to the 
LDA itself and/or numerical accuracy of the total energy. 
With respect to the latter point we note that Fig. || shows 
that the volume-energy curves for the two structures are 
almost parallel near the phase transition. Therefore a 
small change in the total energy for one of the phases is 
expected to cause a large change for the value of Pc; if 
the total energy of one of the two phases is shifted by 
1 mRyd, the value of Pc changes by 20 GPa. 

Figure m shows the pressure dependencies of the atomic 
volume Va, the lattice constants a and c, and the bond 
lengths ri and r2 which are defined as the nearest neigh- 
ber (n.n.) and next n.n.(n.n.n) atomic distances, re- 
spectively. Note that in the bcc phase c 
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and 



ri = r2 = -^a. The filled circles and filled squares 
indicate the results of the present calculations and the 
open squares and opjen inverse triangles represent the ex- 
perimental values.ll3 The solid and the dashed vertical 
lines indicate the boundary of the phase transitions de- 
termined theoretically by us and experimentally by Aka- 
hama et.al!t3, respectively. 

As seen from Fig. ^ the volume variation as a function 
of pressure below 120 GPa (/3-Po) and above 150 GEa, 
(bcc) shows good agreement with the observations .1121 
The obtained pressure dependence of a, c and the bond 
lengths of the (3-Po phase agree well with the experimen- 
tal results.tS Furthermore, the volume reduction at the 
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transition from /3-Po to bcc is estimated to be 0.06 A'^, 
which is in good agreement with the experimental volume 
(about 0.08 A3). Ill 



III. LATTICE DYNAMICS, 
ELECTRON-PHONON INTERACTION AND 
SUPERCONDUCTIVITY 

A. Calculational procedure 

Actual calculational procedures are as follows. We find 
the dynamical matrix as a function of wave vector for a 
set of irreducible q points on a (8, 8, 8) reciprocal lattice 
grid [ 29 points in the IBZ ] for the bcc structure and 
(6, 6, 6) reciprocal lattice grid [ 32 points in the IBZ ] 
for the /3-Po structure. The {I,J,K) reciprocal lattice 
grid is defined in a usual manner: qijk — {i/I)Gi + 
{j / J)G2 + {k/ K)G3, where Gi,G2,Gy, are the primitive 
translations in the reciprocal space. 

The self-consistent calculations are performed for ev- 
ery wave vector with use of the following basis set and 
criteria. We use 3K-spfi-LMT0 basis set (27 orbitals) 
with the one-center expansions performed inside the MT 
spheres up to Emax = 6. In the interstitial region the 
basis functions are expanded in plane waves up to the 
cutoff corresponding to 134 (110), 176 (170), and 320 
(320) plane waves per s, p, and d orbitals for bcc (/3- 
Po) structure, respectively. The induced charge densities 
and the screened potentials are represented inside the 
MT spheres by spherical harmonics up to imax = 6 and 
in the interstitial region by plane waves with the cutoff 
corresponding to the (16, 16, 16) fast-Fourier-transform 
grid in the unit cell of direct space. The fc-space integra- 
tion needed for constructing the induced charge density 
and the dynamical matrix is performed over the (16, 16, 
16) grid [ 145 points in the IBZ ] for the bcc structure 
and (12, 12, 12) grid [ 185 points in the IBZ ] for the /3- 
Po structure, which is twice denser than the grid of the 
phonon wave vectors q. The integration is performed 
also by the improved tetrahedron method. However, the 
integration weights for the k points at these grid have 
been found to take precisely into account the effects aris- 
ing from the Fermi surface and the energy bands. This 
is done with help of the band energies generated by the 
original FPLMTO method at the (32, 32, 32) grid [ 897 
points in the IBZ ] for the bcc structure and (24, 24, 24) 
grid [ 1313 points in the IBZ] for the /3-Po structure. This 
procedure allows us to obtain better convergence results 
with respect to the number of k points. 

For calculation of the electron-phonon coupling the 
corresponding fc-space integrations are more sensitive 
than the dynamical matrices to the number of sampling 
fc-points. It has been performed with the help of the (32, 
32, 32) grid for bcc and (24, 24, 24) for /3-Po in the IBZ 
by means of the tetrahedron method. 



The superconducting transition temperature is cal- 
culated by using Allen-Dynes formula which is derived on 
the basis of the strong coupling theory of phonon mecha- 
nism. Instead of describing the details of the strong cou- 
pling theory, here we give only the necessary equations 
to calculate T^. In the following we completely obey the 
description of the reference.El 

For the electron-phonon spectral d^tribution functions 
a^Fiu)), we employ the expressiont^ in terms of the 
phonon linewidths 7^1^ 
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where iV(eF) is the electronic density of states per atom 
and per spin at the Fermi level. When the energy bands 
around the Fermi level are linear in the range of phonon 
energies, the linewidth is given by the Fermi "golden 
rule" and is written as follows: 

kjj' 

(2) 

where g'^'^qji f,j is the electron-phonon matrix element, 
and conventionally written in the form 
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where kj denotes the one-electron basis 'i'^j and S'^^Vcs 
is the change in the effective potential induced from a 
particular qv phonon mode. Precisely speaking, the 
electron-phonon matrix element must be corrected for 
the incompleteness of the basis functions, but we do not 
discusp-it here. The expression of Tc derived-by AUen- 
Dynesli3 by modifying the McMillan formulall3 is given 
as 
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Usually A is called the dimensionless electron-phonon 
coupling constant, wiog the logarithmic -averaged phonon 
frequency and /i* the effective screened Coulomb repul- 
sion constant whose value is usually taken to be between 
0.1 and 0.15. 

In case of monatomic metals A can be expressed also 
in the following form: 



A = 



N{ev){P) 



V_ 



(7) 



3 



where M is the mass of the atoms and (w^) denotes the 
average of squared phonon frequencies which is given as 
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Further {P) represents the Fermi surface average of 
squared electron-phonon couphng interaction which is 
defined by 
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and 77 = N{eF){P) is called the Hopfield parameter. 



B. Results for bcc Se 

We first calculated the phonon dispersion curve along 
the high symmetry line (FN) for bcc Se at different 4 
volumes (pressures), 10.37 (214.2 GPa), 11.11 
(165.6 GPa), 11.85 (128.6 GPa) and 12.59 A^ 
(102.59 GPa). The resuhs are shown in Fig. |. 

As the pressure decreases, the overall tendency of de- 
crease of phonon frequency is seen. In particular, the 
frequency softening is remarkable for one of the trans- 
verse modes (shown by the solid curve), and this mode 
exhibits a notable phonon anomaly, i.e., a dip in the mid- 
dle-pf the line. The same phonon anomaly is obtained in 
SJl3 This softening of the transverse mode does not cause 
directly the bcc — *■ /3-Po transition with decreasing pres- 
sure because both of the /3-Po and bcc phases haite one 
atom per unit cell. However Zakharov and Cohenllfl have 
pointed out that it plays an important role in changing 
the coordination number from eight to six during the bcc 

/3-Po transition. 

F. Mauri et al. have performed ab initio linear- 
response calculation for lattice dynamics of bcc Te un- 
der pressures.H They reported the same anomaly for the 
transverse mode along the FN line and found that with 
decreasing pressure the phonon frequencies in the mid- 
dle of the FN line become imaginary in a pressure region 
where the /3-Po structure is stable. In our calculation 
complete softening of the transverse mode has not been 
observed even at 100 GPa where the /3-Po structure is 
stable. Complete softening may occure at even lower 
pressures. 

Figure ^ shows the pressure dependence of the phonon 
dispersion along several symmetry lines and the phonon 
density of state (DOS) calculated at three volumes (or 
pressures) 11.85 A^ (128.6 GPa), 11.41 A^ (149.6 GPa) 
and 11.11 A^ (165.6 GPa). It is noted that except along 



the FN line all the phonon frequencies soften linearly with 
decreasing pressure (or increasing volume). 

By using the Allen-Dynes formula we have estimated 
the superconducting transition temperature Tc of bcc Se 
at three pressures: 128 GPa, 150 GPa and 166 GPa. In 
Table || we give the values of calculated Tc together with 
DOS at the Fermi level N{eF), the Hopfield parameter 
77, the logarithmic average frequency wiog, the average of 
squared phonon frequencies (w^) and the electron phonon 
coupling constant A. With decreasing pressure the value 
of wiog decreases while the value of A increases, but the 
rate of change of A exceeds that of wiog. As a result 
the value of Tc increases considerably with decreasing 
pressure. Since A can be expressed by 



A 



the frequency softening (decrease of (w^)) is considered to 
contribute to the increase of A with decreasing pressure. 

In order to obtain a more physical insight into the char- 
acteristic pressure dependence of Tc we consider mode 
and wave-vector dependencies of the phonon linewidths 
7q^ along the symmetry lines. Figure |^ shows that ^q^ 
is almost independent of pressure except for the longitu- 
dinal mode along the FH line and one of the transverse 
modes along the FN. 

With decreasing pressure, 7qj/ of the longitudinal mode 
along the FH decreases whereas that of the transverse 
mode along the FN line increases considerably. Gen- 
erally speaking, a large phonon linewidth increases the 
dimensionless electron-phonon coupling A. Therefore, it 
is expected that the transverse mode along the FN line 
plays an important role in giving rise to the characteristic 
pressure dependence of Tc. 

To clarify the role of the transverse mode along the FN 
line in more detail we have calculated a quantity a'^[uj) 
defined by 



(10) 



where a^F{Lo) is the spectral function and D{lu) denotes 
the phonon density of states. We consider that by in- 
specting the frequency dependence of a^(w) we can dis- 
cern which phonons make dominant contributions to the 
dimensionless electron-phonon coupling A. Fig ^ shows 
the calculated a^(a;) as a function of frequency for three 
pressures. The peak around 2 THz originates from trans- 
verse phonons along the FN line and the peak around 
7^10 THz from longitudinal phonons along the FH line. 
As seen from the figure, both the peaks move towards 
the lower frequency side with decreasing pressure. It 
should be noted, however, that the magnitude of 0^(0;) 
around 2 THz increases remarkably with decreasing pres- 
sure whereas the magnitude of a^luj) around 7'--^10 THz is 
less dependent on pressure. Therefore, we can say again 
that transverse phonons in the middle of the FN line 
make a dominant contribution to A. 
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Combining all of the above results we conclude that the 
origin of remarkable increase of Tc of bcc Se with decreas- 
ing pressure is mainly attributed to the phonon anomaly 
(the remarkable frequency softening) in the middle of the 
FN line. 



C. Results for /3-Po Se 

In this section we calculate the superconducting tran- 
sition temperature of /3-Po Se. To see the pressure de- 
pendence of Tc in /3-Po Se we have calculated Tc at pres- 
sures 103.1 GPa and 118.2 GPa with two sets of lattice 
constants: one is the lattice constants evaluated by cal- 
culation and the other is those determined by experiment 
and given in Table |l|. Figure || shows the electronic dis- 
persion curves and the density of states calculated for 
103.1 GPa. 



The calculated Tc are given in Table III together with 
the electronic DOS at the Fermi level N{eF), the Hopfield 
parameter 77, the logarithmic average frequency LJiog, the 
average of squared phonon frequencies (w^) and the elec- 
tron phonon coupling constants A. Tc strongly depends 
on which sets of lattice constants is used. For the lattice 
constants estimated by calculation' the magnitude of Tc 
is larger and increases considerably with decreasing pres- 
sure. For lattice constants determined by measurements, 
on the other hand, the magnitude of Tc is smaller and 
depends little on pressure. 

The logarithmic average frequencies wjog are larger for 
the experimental lattice constants. The electron-phonon 
coupling A, on the other hand, is larger for the theoret- 
ical lattice constants and furthermore depends consider- 
ably on pressure, which gives higher and more pressure- 
sensitive transition temperatures Tc for the theoretical 
lattice constants. 

In order to clarify the origin of the different magnitude 
and the different pressure dependence of Tc for different 
sets of lattice constants we have calculated the phonon 
density of states Djuj), the spectral function a'^F{u!) and 

{lu) defined by Eq. |lo|. The results are shown in Fig. |l^. 
The magnitude of a^a>) for the theoretical lattice con- 
stants is larger than that for the experimental lattice 
constants throughout almost the whole frequency range. 
In particular the magnitude of a peak in a'^{uj) around 
2 THz obtained for the calculational lattice constants 
is remarkably enhanced compared to that obtained for 
the experimental ones, thus the magnitude of A and Tc 
take large values. For the experimental lattice constants, 
a^(coi) depends little on pressure and hence the magni- 
tude of A and Tc arc also less dependent on pressure. 

Fig. |ll| shows the values of Tc calculated for /3-Po and 
bcc Se as a function of pressure. If we adopt the experi- 
mental lattice constants for /3-Po Se, the superconducting 
transition temperature Tc is almost pressure independent 
in the /3-Po phase and there is a large jump in Tc at the 
transition from /3-Po to bcc. 



IV. CONCLUSION 

We have performed FPLMTO calculation for the /3- 
Po and the bcc structures of selenium, and succeeded 
in reproducing the phase transition from /3-Po to bcc as 
observed by experiments. The obtained pressure depen- 
dencies of lattice parameters agree fairly well with the 
experimental results. The estimated transition pressure 
is 120 GPa. It improves the previous theoretical cal- 
culations, but is still lower than the experimental value 
150 GPa. 

We have investigated the superconductivity of /3-Po 
and bcc Se by calculating the lattice dynamics and 
electron-phonon interaction with use of LR-FPLMTO 
method. For bcc Se we have found that the frequency 
softening is remarkable for one of the transverse modes 
in the middle of the FN-line, and this mode exhibits a 
notable phonon anomaly. The calculated superconduct- 
ing transition temperature Tc increases considerably with 
decreasing pressure, which is mainly attributed to the 
phonon anomaly in the middle of the FN line. In the /3- 
Po structure, on the other hand, Tc is smaller compared 
to that of the bcc phase and is almost independent of 
pressure. Finally we predicte a discontinuous jump in Tc 
at the transition from /3-Po to bcc. 
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FIG. 1. The band structure and the density of states of Se 
at atomic volume y=14.82 A: (a) /3-Po type structure and (b) 
bcc structure. The horizontal line denotes the Fermi level. 
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FIG. 2. The total energy of /3-Po Se as a function of c/a 
for the fixed atomic volume. The inset indicate the results 
for volumes near the phase transition. All the energies are 
referenced to the bcc structure, i.e. AE=E-Ei,cc. 
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FIG. 3. The total energy calculated as a function of vol- 
ume for the /3-Po and the bcc structures of Se. The lines are 
obtained by using Murnanghan's equation of state. 
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FIG. 4. The pressure dependence of atomic volume, lattice 
constants and bond length of selenium. The filed circles and 
filled squares indicate the results of the present calculations 
and the open squaree-and open inverse triangles represent the 
experimental valuesJij The solid and the dashed vertical lines 
indicate the boundary of the phase transition from /3-Po to 
bcc, determined theoretically by us and experimentally by 
Akahama et.al, respectively. 
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FIG. 5. The phonon dispersion for bcc Se. The chain line 
denotes longitudinal mode, and the solid and dashed lines the 
transverse modes. 
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FIG. 6. The phonon dispersion and phonon density of state 
(DOS) for bcc Se. The chain line denotes longitudinal modes. 
The solid and dashed lines denote transverse modes. 
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FIG. 7. The phonon line width calculated for three 
difi'erent pressures. The chain line denotes the results for the 
longitudinal mode, and the solid and dashed ones those for 
the transverse modes. Note that the transverse modes are 
degenerate along the F-H-P-F line. 
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FIG. 8. The frequency dependence of q(cj)^ obtained for 
three pressures. 



0.5 r 



0.0 ■ 



-1.0 ■ 







-J 












V 

























0.5 



0.0 



-0.5 



-1.0 



-1.5 



Z P AT 
Wave vector 



FO 2 4 6 8 10 
DOS [ state/Ryd./atom/spin ] 

FIG. 9. The band structure and the DOS of /?-Po Se at 
atomic volume y=12.59 A (103.1 GPa.). The dashed curves 
denote the results obtained with use of lattice constants es- 
timated by calculation and the solid ones theKe with use of 
lattice constants determined by experiments. 113 The horizon- 
tal line denotes the Fermi level. 
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FIG. 10. The phonon density of states D(uj), the spectral 
function q^F{uj) and a'^{u>) defined by Eq. nO calculated for 
pressures 103.1 GPa (left hand side) and 118.2 GPa (right 
hand side) with use of the theoretical and experimental lat- 
tice constants. The dashed curves represent the resuhs for 
the theoretical lattice constants and the solid curves for the 
experimental lattice constants. 
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FIG. 11. The superconducting transition temperature Tc 
of Se as a function of pressure. The closed circles and closed 
triangles denote the computed values of Tc with /x*=0.12 and 
/x*=Q.10, respectively. For the /3-Po structure the open cir- 
cles and open triangles represent the results for the lattice 
constants estimated by calculation and the closed circles and 
closed triangles those for the lattice constants determined by 
experiments. 



TABLE I. The electronic DOS at the Fermi level iV(£p), 
the Hopfield parameter t], the logarithmic average frequen- 
cies wiog, the average of squared phonon frequencies {lo^), the 
electron phonon coupling constant A and the superconducting 
transition temperatures Tc calculated as a function of pressure 
for bcc Se. The two values for Tc correspond to two different 
values of ^* (0.10 and 0.12). The units of N{eF) and are 
state/Ryd. /atom/spin and Ryd./a^, respectively. 



0.3 


P (GPa)iV(eF) 


V 


c^iog (K) 


{co') (K^) 


A 


Tc (K) 


0.2 


128.6 5.45 


0.20 


224.73 


291.95^ 


0.83 


11.29, 9.90 


0.1 


149.6 5.24 


0.21 


248.03 


316.41^ 


0.73 


9.53, 8.11 


0.0 
1.2 


165.6 5.10 


0.21 


264.62 


335.11^ 


0.66 


8.03, 6.64 



TABLE II. Two sets of lattice constants c/a and a of /3-Po 
Se determined by calculation and experiments. 



[GPa] 


103.1 calc. 


103.1 Exp. 


118.2 calc. 


118.2 Exp. 


c/a 


0.71 


0.75 


0.67 


0.74 


a (A) 


7.471 


7.314 


7.504 


7.255 



TABLE III. The electronic DOS at the Fermi level A''(£f), 
the Hopfield parameter 77, the logarithmic average frequen- 
cies oJiog, the average of squared phonon frequencies {lo'^), 
the electron phonon coupling constant A and the supercon- 
ducting transition temperatures Tc of /3-Po Se calculated for 
103.1 GPa and 118.2 GPa. The upper two lines show the 
results obtained with use of the lattice constants estimated 
by calculation and the lower two lines those with use of the 
lattice constants determined by experiments. The two val- 
ues for Tc correspond to two different values of fi* (0.10 and 
0.12). The units of N{eF) and rj are state/Ryd. /atom/spin 
and Ryd./a^, respectively. 



P (GPa 


)NieF) 


V 


ujiog (K) 




A 


Tc (K) 


103.1 


5.23 


0.19 


192.98 


264.72^ 


0.92 


11.74, 10.49 


118.2 


5.50 


0.19 


204.84 


279.76^ 


0.82 


10.10, 8.84 


103.1 


4.89 


0.16 


250.11 


311. 71"^ 


0.58 


5.14, 4.04 


118.2 


4.85 


0.17 


255.88 


324.58^ 


0.57 


5.01, 3.91 
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